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Abstract. The RBC and UKQCD collaborations have been investigating hadron physics in numer- 
ical lattice quantum chromodynamics (QCD) with (2+1) flavors of dynamical domain wall fermions 
(DWF) quarks that preserves continuum-like chiral and flavor symmetries. The strange quark mass 
is adjusted to physical value via reweighting and degenerate up and down quark masses are set as 
light as possible. In a recent study of nucleon structure we found a strong dependence on pion mass 
and lattice spatial extent in isovector axialvector-current form factors. This is likely the first credible 
evidence for the pion cloud surrounding nucleon. Here we report the status of nucleon structure cal- 
culations with a new (2+l)-flavor dynamical DWF ensembles with much lighter pion mass of 180 
and 250 MeV and a much larger lattice spatial exent of 4.6 fm. A combination of the Iwasaki and 
dislocation-suppressing-determinant-ratio (I+DSDR) gauge action and DWF fermion action allows 
us to generate these ensembles at cutoff of about 1 .4 GeV while keeping the residual breaking of 
chiral symmetry sufficiently small. Nucleon source Gaussian smearing has been optimized. Prelim- 
inary nucleon mass estimates are 0.98 and 1.05 GeV. 
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DISCOVERY OF PION CLOUD? 

In a recent numerical lattice-QCD study [1,2] the RIKEN-BNL-Columbia (RBC) and 
UKQCD collaborations reported a strong dependence on pion mass, m n , and lattice spa- 
tial extent, L, in isovector axialvector-current form factors of nucleon. The dependence is 
most evident in the isovector axial charge, gA, as summarized in Fig. 1. The dependence 
is much stronger than had been observed in non-unitary calculations with DWF quarks 
using either quenched [3] or rooted-staggered fermion [4] ensembles. As the dependence 
appears to scale with a dimensionless quantity, m n L, the product of calculated pion mass, 
m K , and the lattice spatial extent, L, a likely explanation for the dependence is that sig- 
nificant part of the nucleon isovector axialvector current is carried by the expected, but 
never seen, "pion cloud" surrounding the nucleon. If confirmed, this calculation may be 
the first concrete evidence of such a pion cloud. 

Indeed similar strong dependence on pion mass and lattice spatial extent is seen in 
other axialvector-current form factors but not in the conserved vector-current ones [2] or 
in low moments of the structure functions [5], supporting the pion-cloud interpretation. 
However, since the calculations so far have only been carried out at single lattice cut 
off, a~ l , of about 1.7 GeV, two lattice spatial extent, L, of about 1.8 and 2.7 fm, and 
relatively heavy pion mass down to 330 MeV, it seems premature to conclude that the 
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FIGURE 1. Nucleon isovector axial charge, gA, calculated in numerical lattice QCD with DWF quarks 
plotted against a dimensionless variable, m K L, the product of calculate pion mass, m n , and lattice spatial 
extent, L. Results from RBC/UKQCD unitary (2+l)-flavor dynamical quark calculations show scaling in 
this variable with strong dependence. In copmarison, non-unitary DWF calculations with either quenched 
or MILC ensembles show only weak, if at all, dependence. 



dependence is caused by the pion cloud and so we have discovered it. This provides a 
good motivation for the new study being reported here with pion mass as light as 180 
MeV and lattice spatial extent of about 4.6 fm. 



RBC/UKQCD (2+l)-FLAVOR DYNAMICAL DWF LATTICE QCD 

For the past several years the RBC and UKQCD collaborations jointly generated two 
sets of (2+l)-flavor dynamical domain- wall fermions (DWF) lattice QCD ensembles at 
lattice cut off, a -1 , of about 1.7 [6] and 2.2 GeV [7]. We adjust the strange quark mass 
to the physical value by reweighting [8, 7]. The degenerate up and down quark masses 
are set at several values, as light as possible, in the range corresponding to pion mass 
from about 420 MeV to 290 MeV. In contrast to most other lattice fermion methods, 
the DWF formulation provides good control of chiral and flavor symmetries and non- 
perturbative renormalization. This allows us to use meson observables such as pion and 
kaon mass, m K and m#, and decay constants, f n and fx, obtained from these ensembles, 
for a detailed study of meson physics in combined chiral and continuum limit [7]: we 
first use pion, kaon, and Q-baryon mass to set the physical quark mass and lattice scale 
and then obtain predictions for other observables such as pion and kaon decay constants 
with a few % accuracy. 

At this level of high accuracy our predictions is now not limited by our statistics 
but by poor applicability of chiral perturbation or other chiral extrapolation ansatz from 
relatively heavy pion mass we use. Thus the RBC and UKQCD collaborations are now 
jointly generating a third set of (2+l)-flavor dynamical domain- wall fermions (DWF) 
lattice QCD ensembles at lower pion mass of 180 and 250 MeV [8, 9]. This is made 



possible by using a new combination for gauge action of Iwasaki and dislocation- 
suppressing-determinant ratio (DSDR) [10, 11, 12] that allows a lower cut off of about 
1.4 GeV while keeping the residual breaking of chiral symmetry sufficiently small. The 
lower lattice cut off also allows large spatial extent of about 4.6 fm. 

Such a study, however, requires good understanding of chiral and finite-size correc- 
tions to the observables which are unfortunately missing for baryons in general. For this 
reason we are not yet attempting a similar combined chiral and continuum limit study 
for baryons at this moment. Rather, we are refining our fixed-cutoff study of nucleon to 
improve our understanding of its chiral behavior. We have published our findings in nu- 
cleon structure from the 1.7 GeV ensembles in three publications, Refs. [1, 2, 5], which 
we summarize in the following. As we discovered unacceptably large and distorting de- 
pendence on the scaling parameter m n L of axialvector current form factors, especially 
at m % L — 4.5, we did not attempt a similar study for our 2.2 GeV ensembles to save 
our computing resources, and skipped to the new I+DSDR DWF ensembles with lighter 
pion mass and larger lattice spatial extent: the m K L of these ensembles are respectively 
at 4.2 and 5.8 and should allow us to better understand the observed dependence. 



NUCLEON STRUCTURE AT 1.7 GEV 

I here summarize the nucleon structure calculations [1, 2, 5] by RBC and UKQCD 
collborations using the 1.7-GeV dynamical (2+l)-DWF ensembles [6]. 

In an earlier study with two dynamical DWF flavors [13], we identified an important 
source of systematic error in lattice-QCD numerical calculation of nucleon structure, 
namely excited-state contamination. As the amount of contamination varies depending 
on the shape of source smearing, we need to optimize the combination of source smear- 
ing and source-sink separation in order to filter out the excited-state contamination while 
maintaining reasonable statistical signal. If we choose too long a separation then even 
the ground state decays and no signal is obtained. For our choice of Gaussian source 
smearing width of 7 lattice units, a source-sink separation of about 12 was optimal. 

After making this important adjustment we discovered the isovector axial charge 
strongly depends on pion mass and lattice size, as was discussed earlier with Fig. 1. 
The dependence seems to scale in a single parameter, the product x = m K L of the pion 
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FIGURE 2. In 2-flavor dynamical DWF study [13], systematic difference was identified arising from 
excited-state contamination. Source smearing and source-sink separation must be optimized to filter out 
the excited state while maintaining reasonable statistical signal. 
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FIGURE 3. Isovector vector current form factors. Left: Dirac form factor mean square radius. Center: 
Pauli form factor mean square radius. Right: isovector anomalous magnetic moment. The former two show 
no sign of expected singularity at m K = 0, and significantly undershoot experiments. The latter seems in 
rough agreement with the experiment. 



mass m K and lattice size L. Though our dynamic range in the scaling parameter is 
rather narrow to distinguish various ansatze on functional form, f(x), of this scaling, 
and so cannot yet clarify if this is indeed the pion cloud, fitting to various forms such as 
(inverse volume) or e~ x / ' y/x (pion cloud) is possible, and results in an estimate 
of gA = l-20(6)stat(4) S yst- This strong dependence on m K L is observed also in other 
axialvector-current form factors. In order to drive the systematic error arising from this 
dependence below 1 %, we would need m K L of 6 to 8, or L of about 5 fm for m n of 300 
MeV and 10 fm for physical pion. 

In contrast to the axialvector-current form factors, the vector-current ones do not show 
dependence on the lattice size even at the lightest pion mass of about 330 MeV. Mean 
squared radii of Dirac and Pauli form factors are obtained (see Fig. 3) without chirally 
expected singular behavior in m|, and undershoot the experiments. On the other hand 
the anomalous magnetic moment is in rough agreement with the experiment (Fig. 3). 
These observations are confirmed by a LHP study using our 2.2-GeV ensembles [14]. 

In low moments of isovector structure functions, such as quark momentum fraction, 
(x)u-d> or helicity fraction, (x)Au-Ad, no dependence on m K L is seen either. In a naturally 
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FIGURE 4. Naturally renormalized ratio, (x) u -dj (x)au-Ac1, of the isovector quark momentum and 
helicity fractions, does not show dependence on mass nor lattice size and agree with experiment. 
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FIGURE 5. Isovector quark momentum fraction, (x) u -d (left), and helicity fraction, (x)au-M (right). 
Non-perturbatively renormzlized and run to MS 2 GeV, z MS ( 2GeV) = 1.15(3), and 1.15(4), respectively. 
Neither shows volume dependence between larger, L — 2.7 fm, and smaller, 1.8 fm, lattices. Both show 
interesting trend down to the experiment at the lightest mass. 



renormalized ratio, (x) M _^/ (x)Au-Ad> of the two fractions, no dependence on the mass, 
m q oc m\, is seen either and the agreement with experiment is excellent (see Fig. 4). 
No dependence on L is seen in respective fractions either (Fig. 5). Both show interesting 
trending down toward the experiment, motivating us to calculate these quantities at lower 
pion mass. 



STATUS AT 1.4 GEV 

Our studies of (2+l)-flavor dynamical DWF lattice QCD in (~ 3fm) 3 spatial box and 
down to about 300 MeV pion mass [7] clearly point to the need of numerical calculations 
at yet lighter pion mass, both for combined chiral and continuum limit study in the 
meson sector and for better understanding of chiral behavior in the baryon sector. Thus 
RBC and UKQCD collaborations started to generate a new set of ensembles with pion 
mass about 180 and 250 MeV [9]. As the lower pion mass demands larger lattice spatial 
extent, the new ensembles are generated with about 4.6 fm spatial extent that translates 
to the scaling parameter m n L of above 4 for the lighter mass and almost 6 for the heavier. 
This is made possible by a newly developed gauge action with an appropriate addition 
of dislocation-suppressing-determinant-ratio (DSDR) term [10, 11, 12]. We have about 
1900 hybrid Monte Carlo time units for the heavier ensemble and 1200 for the lighter, of 
which first 600 and 500 respectively are discarded for thermalization. We analyze every 
eight time unit with four evenly separated nucleon sources each. 

We use Gaussian smearing [15, 16] for nucleon source to optimize the overlap with 
the ground state and compared the cases for widths 4 and 6 lattice units. We found width 
6 is better for both pion mass [17]. We quote preliminary nucleon mass estimates of 
0.721(13) and 0.763(10) lattice units which correspond to about 0.98 and 1 .05 GeV with 
another preliminary estimate for the lattice cut off of 1.368(7) GeV. As are summarized 
in Fig. 6, these ensembles are filling the gap toward the physical point. We look forward 
to reporting nucleon form factors and low moments of structure functions soon. 
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FIGURE 6. Nucleon mass from the RBC and UKQCD (2+ 1 )-flavor dynamical DWF ensembles plotted 
against corresponding pion mass squared. 124 and 132 are from ensembles with simple Iwasaki gauge 
action while ID32 are with new Iwasaki+DSDR gauge action. 
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